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ABSTRACT 


This paper describes a procedure for experimentally determining, in 
terms of the particle motions, tht shapes of the low order acoustic modes 
in enclosures. The procedure is based on finding differentiable functions 
which approximate the shape functions of the low order acoustic modes when 
these modes are defined in terms of the acoustic pressure. The differentiable 
approximating functions are formed from polynomials which are fitted by a 
least squares procedure to experimentally determined values which define the 
shapes of the low order acoustic modes in terms of the acoustic pressure. 

These experimentally determined values are found by a conventional technique 
in which the transfer functions, v;hich relate the acoustic pressures at an 
array of points in the enclosure to the volume velocity of a fixed point 
source, are measured. The gradient of the function which approximates the 
shape of a particular mode in terms of the acoustic pressure is evaluated to 
give the mode shape in ‘terms of the particle motion* The procedure was 
tested by using it to experimentally dete!<nnine, in terms of the pai ticle 
motions, the shapes of the low order acoustic modes in a small rectangular 
enclosure. Details of the experimental technique used and some of the 
results obtained' are given. 
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INTRODUCTION 


r 


The excitation of one or more of the low order acoustic modes in an 
enclosure can generate in the enclosure sound pressures with distinct tonal 
qualities which can annoy the occupants of the enclosure.. A commonly 
experienced example of this situation is the "cavity boom” which sometimes 
occurs in the cabins of passenger cars. Cavity boom arises when a low 
order acoustic mode, which may be coupled to a complaint portion of the 
enclosing structure, is excited. Usually the cavity boom is excited by the 
engine and it often occurs v/heri the engine firing frequency coincides with 
the frequency of a low order mode. Many studies, both theoretical and 
experiiaental, have been made to investigate the low order acoustic modes 
in passenger Annoying sound presstires associated with the presence 

of low order acoustic modes arise in the cabins of agricultural and industrial 
machines!^. A similar effect Ccin occur in the cabins of ships and aircraft. 

It is usually a simple task to experimentally determine the frequencies ^ 
of the. low order acoustic modes of an enclosure. However/ it is not a simple, 
task to experimentally determine the mode shapes, particularly if they vary 
with more than or#^ spatial co-ordinate and if the modeshapes are required 
in terras of the particle motions. A knowledge of the modeshapes, eind in 
particular the modeshapes in terras of the particle motions, is useful in 
devising appropriate noise control measures. 

Over the last decade there have been considerable developments in 
experimental raethodls for determining the modeshapes of structures®'®. The 
ability to describe the motion of a structure when it is vibrating in a 
particular mode is one of the most useful developments which has occurred, 

The ability to describe the motion, has been facilitated by the fact that 
a motion, transducer, usually an accelerometer, is one of the two transducers 



which aice used to generate the electrical signals from which the transfer 
functions, used in the experimental modal analysis of the structure, are 
obtained. Although considerable developments have been made in experimental 
methods for determining the modeshapes of structures,, it is only in recent 
years that attempts have been made to experimentally determine the shapes of 
acoustic modes* One possible reason for the relative lack of progress 
in the acoustic case is that the most conveniently measured acoustic 
queintity, the acoustic pressure, ceinnot be used directly to df/termine the 
particle motion associated with a particular mode. 


This paper describes a procedure whereby the particle 'motions 
associated with* lov; order acoustic modes in enclosures can be derived from 
tj^ansfer function measurements which relate tho acoustic pressures at an 
array of points in an enclosure to the volume velocity of a fixed point 
source in the enclosure. The theory of the procedure is described and than 
the results of a numerical evaluation of the procedure are given. Finally, 
a description of how the procedure was used to experimentally determine the 
low order acoustic modes in a small rectangular enclosure is given. 






I. DESCRIPTION OF THE METHOD 


OF pma ^ 


Let an 3«yz axes system be located in the enclosiire as shovm in Pig, 1. 
The axes are not necessarily orthogonal, suppose that LLX MMX NN transfer 
functions which relate the acoustic pressures at the ILX mmx nn points 
shown in Pig. 1 to the volume velocity of a fixed point source are measured. 
The points shown in Fig. 1 are on lines par’allel to the X/ y and z axes. 
The distances between the points on a line parallel to a particulcu: axis need 
not be equal but the distances between corresponding points on different 
lines are equal. Suppose, that by a conver.tional modal analysis method, 
such as the method described by Nieter and Singh®, the acoustic pressure at 
the point Xxi , yjj , z^k unit volume velocity of the source is measured 
as Pr (x:ii,yjj,zjyQ) in the r^^ mode whose angular frequency is Wr • This 
measured quantity is composed of the true value, denoted Pr* (^tix/Yjj/ZKK) 
and an error term, denoted . Thus: 


Pr^^lI^yjJ'ZKK) “Pr*<i'IlfyjJ/2KK> + E 


The IjLX MMX nn values of derived from the LLX mmx nn 

measured transfer functions cein be used to determine a function which- 

approximates Pr* (x,y,z) . Let this approximating function, which is denoted 
* 

Pj-(x,y,z) , be given by: 

Pj.{'x,y,z) = cCl+XjX-t... +Xjx^”^ +... +XxX^ X 
(1 -HYgy + . . . +Yjy‘^“^ + . . . +yMy”~^ > ^ 

(1 -tZjZ + ... +Zj^2 +... +Zj^z (2) 


The constants C , . . . , Xj,... , Yj,... , which define the approximating 

function can be derived from the measured values, Pr , by the 
following procedure. A total of MMX NN polynomials of the form: 


y f + y * X + 

XyJJ^KK 2,JJ,KK 




+ X • X 

ri,jj,KK 



OF POOR QUAU'fV 


can be fitted by a least squares procedure to the LL values of Pr 
along each of the MMX NN linos which are parallel to the x axis . The 
coefficients derived from a least squares fit along a particular 

line can be nor-aalized by dividing all of the coefficients by one of the 
coefficients. If the first coefficient is used as the normalizing coefficient, 
a polyiromial of the following form is obtained: 


1 +. 






Xi,JJ,KK is given by . Ideally, corresi'onding normalized 

0 

coefficients for each of the MMX nn polynomials will be equal. However, since* 
these coefficients are derived from experimental data they will differ and so 
Xj can be estimated by averaging the Xi^jj^lO^ coefficients. Thus Xj is 
given by: 


' MM NN 

Xj = S 2 Xx^jj^Kk/MMXNN 
JJ«X KK«1 


Along some lines parallel to the x axis the experimentally derived values 

small and so the coefficients and hence the noinnalized 
coefficients may not be representative and it is prudent to discard such 
coefficients in the averaging. This can be done by evaluating the **size” of 
each of the fitted polynomials in terms of a measure such as the root mearc 
square value of the fitted polyriomial and then discarding coefficients 
associated with polynondals whose ‘’size^* is, for example, less than the 
average size, • 

A similar procedure involving the fitting of polynomials in the y and z 
directions can be used to determine the normalized coefficients **.yj ...Zk ^ 

Once the normalized coefficients • . • • have been found, the 

coefficient C which appears in Eqn. (2) can b jmd to minimise the sum 
of the squares of the differences between Pr and Prt^IIfYjjrZKK) • 


’i ; 








•OF POO.? Qumr^ 


It is found, in the development of the acoustic wave equation, that 
the acceleration of a gas particle in a particular direction is given by 
the negative of the gradient of the pressure in that direction divided by 
the density of the gas. Thus, for the mode, the amplitude of the particle 

acceleration in the x direction* per un-*.t volttme velocity of the source 
arx(X/y#z!) , can be approximated by; 

. arx(.=,y,z) (3, 

It can be seen frorii this equation that the first partial derivative of 
the approximating function with respect to x must be found. This can be 
readily done v/hen is defined in the form of Eqn. (2). The 

amplitudes of the particle t velocity and the particle displacement can be 
readily found from the cuuplitude of the acceleration. Similarly/ the 
amplitudes of the acceleration/ velocity and displacement in the y and z 
directions can bo estimated from the partial derivatives of the approximating 
function with respect to y and z. 


II. CALCULATIONS WITH THE FITTING PROCEDURE 

» » 

< , 

It is evident that the approximating function Pr(5s.y>z) must approximate 

sufficiently well so that the first partial differentials of 
Pr(x,y,z) and P;u*(x,y,z) are in satisfactory agreement. An indication of 
how many terms should be incorporated in the approximating function can be 
obtained by using a suitable differentiable function to generate data at the 
points Xjj,yjj,Zj^; 11= 1,LI,; JJ= 1,MM; KK= 1,NN , fitting approximating 
functions with various nvimbers of terms and then comparing the first partial 
differentials of the original and the approximating functions. 

A suitable differentiable fxmction is given by; . 
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ORICffM! 

OF POOR QUALITY 


(Attx) COS (roiTy) COS (nirz) (4) 

** 0 jf X f 2 f • • • 

This function, in addition to being yei^dily differentiable, defines 
in terms of the acoustic pressure, the shapes of the modes in a rectangular 
enclosure whose sides have unit lengths in the x, y and z directions. 

There are H>, m and n half waves in the x , y and z directions. 

Fig. 2(a) shows the values of the original and the approximating 

functions along the x axis after an approximating function with 3j»MwN«3 

was fitted to Lba MM» NN» 4 values derived from the original function, Eqn. (4) 

with m» n= 1 . The gradients along this line in the x , y and z directions 

are shown in Pig. 2(b) for the original and the approximating fuiiction*!. Figs. 

* » 

3(a) and 3(b) show the corresponding results for the case of L« M» N» 4, 
lx =MM »NH =5 and )L«=m«n»l . Compeirison of Pigs. 2(b) and 3(b) indicates- 
that if a mode in an enclosure involves a half wave in a given direction, 
the corresponding approximating polynomial should incorporate at least four '• 
coefficients. It is of interest that, for a rectangular enclosure, seven 
modes can occur v;hen there is no more than one half wave in each of the x , y 
cuid z directions. 






Fig. 4(a) shows the values of the original and approximating functions 
along the x axis after an approximating function with L=M«N='6 was fitted 
to LI»» MM= NNi= 7 values derived from the original function, Eqn. (4), with 
S. «m =n =2 . The gra(kients along this line in the x , y and z directions are 
shown in Pig. 4(b); for the original and the approximating functions. Pigs. 

• 5(a) and 5(b) show the corresponding results for the case of L=M=N=7, 

LL =MM =NN =8 and £. =m =n = 2. Comparison of Figs. 4(b) and 5(b) indicates 
tfeat if a mode in an enclosure involves two half waves in a given direction, 
the corresponding approximating polynomial should incorporate at least seven 
coefficients. 
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It lo of lntaro 0 tf that/ for a rectangular onclosurs/ fifteen modes can 


occur when there are no more than two half waves in each of the x , y and z 


directions. 


Although it is possible to fit polynomials to modes which involve more 


than two half v/aves in each direction/ the number of coefficients required in 


the approximating polynomial becomes large as does the number of points at 


which measurements must be made. The number of data points at which 


measurements must be made depends upon the magnitude of the error term B in 


Eqn. (1) relative to the magnitude of the true value of the function. If 


more than two half waves occur in a given direction/ the approximating 


function can be fully or partially constructed from harmonic functions instead 


of the polynomials used in- Eqn. (2}. 


in. EXPERIMENTS 


The previously described procedure was applied to experimentally 


determining the shapes of the particle velocity modes in a small rectangular 


enclosure. These modeshapes can be compared with the theoretically determined 


modeshapes. 


The internal dimensions of the rectangular enclosure and the orientation 


of the xyz axes system are shown in Pig. 6. The acoustic field in the 


enclosure was excited by a horn driver which was located in the x => 0 # y = 0 / 


z =h^ corner of the enclosure. The horn driver, which was driven by noise 


which was white over the 50 Hz to 500Hz range contained a microphone in the 


magnet cavity so that the volume velocity of the source could be measured. 


Details of the source and how it was calibrated are given in Appendix A. 
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Transfer funcUions which related tlie acoustic pressures at points along 
lines parallel to the k # y -®nd z axis to the volume velocity of the source 
were measured. The points, which were equally spaced along the lines, were 
defined by the values Hi* 11, MM ■«7 and NN *9 . The values of Pr 
were determined for the first ton resonance frequencies by the coincident- 
quadrature procedure® . An approximating function was fitted to each of the 
ten sets of data by the previously described procedure. First, all of the 
sets wore fitted with polynomials defined by L» M« N» 6 so that the general 
shape- of the mode in the direction of each axis could bo estimated. The 
appropriate values of the polynomials in each direction for each set of data 
could then be fixed. If there were two half v/aves in a given direction tlie 
polynomial index L, M or N for that direction was set to 7 on the basis of 
the results given in the previous section. If there was one half wave in 
a given direction, the polynomial index L, M or * N for that direction was 
set to 4. If the acoustic pressure amplitude was nominally constant in a 
given direction, the polynomial index L, M or N for that direction was set 
to 2. , 

A sample of the results derived is presented in Figs. 7 and 8. The 
results given in Fig. 7 relate to the first mode, whose frecjuency was 
123.8 Hz. It was found, with reference to Pig. 6, that there was one half 
wave in the x direction and no variaition in the y and z directions. The 
amplitude of the acoustic pressure per unit volume velocity of the source 
varied along the x axis as shown in Fig. 7 (a) . Values derived from the 
approximating function and from the measurements along a line close to the 
axis are shown. The particle velocities derived from the approximating 
function varied along the x axis as shown in Fig. 7(b)- The results given 
in Fig. 8 relate" to the tenth mode whose frequency was 302.8 Hz. It was 
found, with, reference to Fig. 6, that for this mode, there were two half 
( waves in the x direction, one in the y direction and there was no variation 
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In tho z direction. The ninth mode, which Involved two half waves in 
tho z direction and no variation in the x and y directions occurred at 
a frequency of 296.2 llz, only 2% lower than tho frequency of the tenth mode. 
The anplitudes of the acoustic pressure per unit volume velocity of the 

# 

source varied along the x and y axes as shown in Pigs. 8(a) and 8(b). 
Values derived from the approximating functl6n and from tho measuremontn 
along lines close to those axes are shown, The particle velocities derived 
frem tho approximating function ^u:e shown in Pigs. 8(c) and 8(d). 

Although the results given in Pigs . 7 (b ) , 8 (c) and 8 (d) contain the 
expected features of themodeshapes in terms of the particle motions, they 
do nob exhibit the expected results at the boundaries. Thus, for example, 
in Pig. 7 (b) , it can bo seen that although tho particle velocity in the x 
direction is zero at x« 0, it is not zero at x« 1.482m. This feature was 
found to present to a lesser extent in the corresponding results for 

ot^iiijr This situation is probably due to a number of factors. such 

as tho non-rigid boundary of tho enclosure used in the experiments and 
the irregular data which defined the shapes of the modes in terms of the 
acoustic pressure. A comparison of Pigs. 7(a) and 8(a) shows that the 
experimental data for the first mode was in geniiral more irregulcir than that 
for the tenth mode. However, if it is certain that a pcirticular boundary is 
rigid, the least squ^es fitting procedure can be readily modified to make 
the gradient of the approximating function at and normal to that boundary 
zero. The particle velocity at and normal to that boundary v/ill then be 
zero. 
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IV. CONCLUDING REMARKS 

Tha ejtperimental results derived from the procedure described in this 
paper suggest that the procedure can be developed to allow the peirticle 
motions associated with low order acoustic modes in enclosures to be found. . 

The possibility of displaying these modes in animated form as is presently 
donu in displaying the vibrational modes of a structure then exists. 
Undoubtedly/ superior approximating functions than the one defined by Eqn. 

(2) exist eind. several types of approximating functions are currently being 
investigated . . ’ 

^ 'v 

Particular attention is i>eing paid to tho ability of these approximating 
furictions to extract useful information from irregular data. However, the 
approximating funct*ion defined by Eqn. (2) can be rapidly determined with 
basic computafci^^nal facilities, a feature not possessed by other approximating 
functions which v;ere considered. 
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APPENDIX A 


THE CALIBRATED VOLUME VELOCITY SOURCE 


The calibrated volume velocity source was constructed from the driver 

of a horn type loudspeaker. The horn driver was modified- by grinding slits 

into the magnet so that the volume between the diaphragm and the magnet was 

freely connected to the volume around the magnet. Pressure changes in the 

volume around the magnet due to the motion of the diaphragm of the driver 
1 

could be measured by a microphone which was inserted into the volume around 
the magnet. A cross-section of the arrangement is shown in Fig. A.l. Thus 
by measuring the pressure with the microphone, the volume velocity produced 
by the source could be determined. 

It can be easily shown that^ if the total volume enclosed by the 
diaphragm, Vq , changes by a small amount, 5v, the absolute pressure, Pq , 
in this volume chcmges by a small amount, 6P =-6v ^ (Y^o/^o) / in 

the volume behaves in an adiabatic manner. The specific heats ratio is given 
by y . Thus the volume velocity b <■:::) of the source is related to the 
acoustic pressure p(t) in the volume behind the diaphragm by Eqn. (A.l): 


U(t) ^ -(Vo/YPq) ^ dp(t)/dt 


(A.l) 


The volxme enclosed by the diaphragm, Vq , being of very irregular 
shape,, could not be accurately determined and so the calibration was made 
by measuring the acoustic pressure in a small accurately known volume which 
was attached to the throat of the driver along with the voltage generated by 
the microphone which was used to measure the acoustic pressure in the volume 
enclosed by the diaphragm. The volume velocity was derived from the acoustic 
pressure measured in the small accurately known volume by use of an equation 
of the form of Eqn. (A.l). The amplitude of the frequency response function 
relating the volume velocity produced by the source and the voltage generated 




by thn microphone is shown in Pig, A. 2, This frequency response function 
was derived from a transfer function measurement with a white noise 
electrical signal being applied to the coil of the driver. 



CAPTIONS 


Pig. 1 Nomenclaturo for measurement points. , 


Pig. 2 (a) Magnitudes of the original and the approximating ( ) 

functions along the x axis. 

(b) Gradients in the x, y and z directions of the original ( 

and the approximating C — ~) fui ctions along the x axis. 



3(a) 

(b) 

4(a) 

(b) 

5(a) 

(b) 


L “M =»N =3 ; LL »MM «NN =4 ; & =m =n =» 1 . 

* 4 

Magnitudes of the original ( ^) and the approximating ( )’ 

functions along the x cixis. 

Gradients in the x , y and z directions of the original ( ) 

and the approximating (- -) functions along the x axisl 

L =M =N =4 ;LI, =MM =NN =5 ; A =ra =n =1 . . 

Magnitudes of the original ( ) and the approximating ( •') 

functions along the x cixis. 

Gradients in the x, y and z directions of the original ( ) 

and the approximatim ( — ) functions along the x cixis. 

L »M =N ?=6 ; LL =MM =NN =7 ; Z =m =n =2 . 

Magnitudes of the original ( ) and the approximating ( ) 

functions along the x axis. 

Gradients in the x / y and z directions of the original ( ) 

and the approximating ( — ) functions along the x axis. 

L =M =N =7 ; LL =MM =NN =8 ; i =m =n =2 . 
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Pig. 2<a) Magnitudes of the original ( ) and the approximating ( 

functions along the x axis. 

• (b) Gradients in the x , y and z directions of the original (■ 

and the approximating C — — ) functions along the x axis. 
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Pig, 3(a) Macjnitu«3os of tho original zmd the approximating ( 

functions along the x axis. 

(b) Gradients in the x, y and z directions of tho original (• 
and the approximating functions along the x axis. 

. Th “M »4 ; LI# "MM "IW “5 ; f "m »*n »1 . . 









Pig. 3 (a) Magnitudos pZ the original (- 
functions 'along tho x axis* 


-) and tho approximating ( 


Gradients in tho x , y and z directions of the original ( 


and the approxitnating (• 


•) functions along tho x Eucis. 
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4 (a) Magnltudoo oC tho original (———*) ami the appi^oyimating 

« * 

functions along the x axis, , * ' 


(b) Gradients in tho x, y and z directions of tho orlginaX ( •— ■) 

and the approximating , functions along tho sf axis. 
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I, «M -N -6 ? tL -MM -NN «7 ; % »>rt -2 . , 
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Fig. 5{«) Magnitudes of the original ( — ) and the approximating ( 
functions along tha x axis, 

(b) Gradients in tho x, y and z direetlons of tljo original (> 
and tho approximating ) functions along tho x axis. 

It »M «N »7 ; bL »MM “NN "8 ; ^ «m «»n “2, . 
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Jf'ig. 7(a) Pressure amplitude por unit volume velocity ^lmplitudo along 

the X axis frcm the measured values (•) and from the approximating 

* . • » » , . • 
function { : — ■) . 


(b) Particle velocity amplitudes pec wit 'volume velocity amplitude 
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Fig, 7(ei> Prossuro amplitude per unit volume velocity amplitude along 
» • 

the X axis fron the measured values {•) and from the approximating 

function ( — ) . . . 

* » » * 

(b) Particle velocity amplitudes per unit 'volume velocity amplitude 
. along the x axis. ■ ' ‘ , 

L =4 ,M BiJ “2 j IX. bIL ,MM =7 #NN =9 , 
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Pig. 8(a) Pressure amplitude per unit volume velocity 
amplitude along the x axis from the measured 
values ( . ) and from the approximating 
function ( ) . * 
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Pig. 8(b) Pressure amplitude per unit volume velocity 

amplitude along '.the y. axis from the measured values 
(.) and from the approximating' function ( ). 


O 

VO 








OF POOR QUALITY 


Pig. 8(c) 


Partxcle velocity' amplitudes per unit volume 
velocity amplitude along the x'axis. 
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Pig, 8 (d) ^ Particle velocity amplitudes per unit volume 

velocity amplitude along the y axis. 
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